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CA.  J.R.  Maxwell  of  the  Environmental  Research  Institute  of  Michigan 
provided  assistance  in  using  and  interpreting  the  tapes. 
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SECTION  I . INTRODUCTION 


Infrared  threat  warning  receivers  offer  the  attractive  possibility 
of  improving  aircraft  survivability  through  sensing  and  detection  of 
emitted  and/or  reflected  radiative  threat  characteristics.  This  is 
especially  true  where  the  threat  uses  passive  sensing  and  homing 
techniques.  Even  in  the  absence  of  active  energy  from  the  threat,  the 
IR  warning  receiver  can  detect  thermal  radiation  from  the  rocket  plume 
or  from  the  elevated  temperature  of  a threat  vehicle. 

The  major  problem  encountered  in  the  development  of  an  IR  warning 
receiver  is  due  to  the  competing  requirements  for  high  sensitivity  and 
low  false  alarm  rate.  Thus  infrared  warning  receivers  have  had  limited 
utilization  due  to  an  inability  to  discriminate  against  certain  back- 
grounds which  produce  false  alarms.  Current  and  past  methods  for  solving 
this  problem  have  required  the  use  of  more  and  more  sophisticated  dis- 
crimination techniques  which  rely  on  spectral,  spatial  and  temporal  char- 
acteristic differences  between  the  threat  and  its  background.  Backgrounds 
which  produce  high  false  alarm  rates  have  been  identified  and  programs 
are  underway  to  provide  high  resolution  spectral  measurements  and  models 
for  more  effectively  discriminating  against  these  background  sources. 

Little  is  known  about  the  probability  of  encountering  the  identified 
classes  of  high  radiance  background  sources.  Thus  it  was  the  objective 
of  this  effort  to  develop  statistical  descriptors  for  background  data 
using  available  measurement  data.  The  descriptors  were  to  show  inter- 
dependence of  radiance  and  area  for  the  background  data. 
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Much  of  the  existing  background  data  has  been  gathered  in  recon- 
naissance efforts  utilizing  banded  radiometric  scanners;  thus,  no  high 
resolution  spectral  data  were  available  for  analysis  at  this  time.  The 
inclusion  of  spectral  dependence  offers  another  dimension  for  comparison 
and  discrimination  and  could  be  incorporated  at  a later  time.  Correlation 
between  spectral  bands  for  existing  data  could  be  accomplished  and  has 
been  investigated  by  others.  [1]  Thus  the  major  emphasis  for  this 
effort  was  on  amplitude  and  area  statistics  for  individual  records  in  a 
given  spectral  band. 

The  following  specific  tasks  are  listed  in  the  statement  of  work. 

1.1  Develop  a basis  for  characterization  of  background  sources  such 
as  stack  burn-off,  stack  emissions,  glint  off  water,  etc.,  which 
will  provide  their  radiant  emittance  and  solar  reflectance  pro- 
files in  the  4.3  micrometer  (C02  band)  spectral  region. 

1.2  Reduce  AFAL  furnished  multispectral  data  to: 

1.2.1  Identify  and  classify  specific  background  sources  with 

4.3  micrometer  spectral  signatures  which  are  indicative 
of  source  temperatures  well  above  the  ambient. 

1.2. 1.1  For  each  class  of  these  sources  determine  the 
distribution  of  the  characterizing  parameter, 
e.g.,  radiant  emittance,  spatial  profile,  etc. 

1.2. 1.2  Develop  statistics  describing  the  composition 
and  distribution  of  the  background  sources  in 
various  scenarios. 

1.2.2  Develop  iso-radiant  emittance  terrain  models  (similar  to 
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SECTION  II.  THEORETICAL  DEVELOPMENT  FOR  STATISTICAL 
DESCRIPTION  OF  BACKGROUND  RADIANCES 

The  essential  features  for  modeling  background  radiances  measured 
by  a sensor  and  relating  that  radiance  to  false  target  probability  are 
shown  in  Figure  1. 

A statistical  model  for  evaluating  false  detection  probability  for 
backgrounds  is  best  developed  in  terms  of  conditional  probabilities.  The 
single  glimpse  probability  of  detection  for  a false  target  is  given  by 

p = p,.p2-p3  n) 

where  P = single  glimpse  detection  probability  for  a false  target 
P-|  = probability  of  detection  given  that  the  false  target  is 
in  the  sensor  field-of-view 

= probability  the  false  target  is  in  the  sensor  field-of- 
view  given  the  false  target  is  within  the  scene 
P3  = probability  the  false  target  is  within  the  scene. 

P3  may  be  further  developed  as  a conditional  probability  in  terms  of 
general  scene  categories  such  as  urban,  suburban,  industrial,  forest, 
desert,  etc. 

P-j  depends  primarily  on  sensor  NEFD  and  other  performance  measures, 
false  target  irradiance  levels  at  the  sensor,  and  detection  threshold. 
This  probability  can  be  expressed  in  concise  terms  showing  dependence 
on  parameters  from  the  target,  the  optical  path  and  the  sensor. 

?2  depends  on  the  sensor  instantaneous  field-of-view,  the  total 
search  angle,  scan  geometry  and  on  the  sensor  relative  motion.  Models 
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topographical  representations)  indicating  the  distri- 
bution and  mix  of  the  different  classes  of  background 
sources. 

1.2. 2.1  Develop  appropriate  algorithms  and  determine 
the  probability  of  encountering  these  back- 
ground sources  as  a function  of  scenario. 

Specific  items  in  the  statement  of  work  were  addressed  according  to 
characteristics  of  available  data  tapes  and  supporting  information.  The 
available  data,  which  are  described  in  Section  3.0,  were  limited  to 
relatively  benign  backgrounds  for  the  infrared  warning  receiver.  Additionally, 
the  tapes  were  received  only  two  months  prior  to  the  contract  completion 
date.  Thus  limited  analyses  were  accomplished,  and  the  emphasis  was  shifted 
to  demonstrating  specific  analysis  methods  which  could  be  applied  to  back- 
ground data  impacting  infrared  warning  receivers  as  it  becomes  available. 

Section  2.0  gives  a brief  theoretical  development  describing  the 
effects  of  various  sensing  parameters  on  a statistical  model  for  background 
radiances.  Section  4.0  describes  the  analysis  algorithms  used  for  showing 
radiance-area  characteristics.  Again,  these  do  not  represent  all  analyses 
which  could  be  implemented;  however,  they  do  provide  useful  descriptors. 

Section  5.0  gives  examples  of  the  data  analysis  results  for  the 
limited  data  that  was  available.  And  finally,  Section  6.0  gives  a summary 
and  conclusions  about  the  implemented  analysis  methods. 

Complete  software  listings  are  given  in  the  Appendix. 
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Figure  1 . Modeling  Features 
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for  P2  will  typically  be  based  on  specific  scan  geometries 
for  a given  f ield-of-view.  Two  extreme  models  are  1)  very 
wide  f ield-of-view  sensor  with  slow  angular  scan  and  2)  narrow 
f ield-of-view,  high  spatial  resolution  with  fast  angular  scan. 
For  the  first  case  P2  will  be  large;  however,  P^  will  also  be 
affected  since  sensor  performance  depends  on  both  the  sensor 
field-of-view  and  the  scan  rate.  For  the  second  case  with 
high  relative  aircraft  movement  there  may  be  regions  within  the 
background  which  are  not  scanned  at  all.  This  possibility 
must  be  considered. 

P^  depends  on  the  general  scene  category,  relative  sun 
angles  for  false  targets  due  to  sun  glint,  and  obscuration 
by  low  radiance  objects  (important  primarily  for  sensing 
at  small  elevation  angles) . A general  development  for  P^ 
is  impossible  because  of  the  many  unknown  parameters  involved. 
The  approach  to  P^  evaluation  will  be  based  on  empirical  data. 

The  above  probabilities  are  defined  for  a single  scan  of 
the  search  volume  and  based  on  irradiance  threshold.  They 
do  not  include  the  effects  of  spatial,  temporal  or  spectral 
discrimination  methods.  Thus  P as  defined  in  (1)  represents 
a worst  case  false  target  detection  probability. 

This  report  examines  P^  by  developing  statistics  for 
empirical,  multi-spectral  sensor  data.  The  following 
sections  describe  the  data,  the  statistics  computed  and 
results  of  statistical  characterization  of  the  data. 


SECTION  III.  DESCRIPTION  OF  DATA  TAPES  TO  BE  ANALYZED 


A single  data  tape  containing  nine  files  (scenes)  of  scanned,  banded 
radiometric  data  in  three  spectral  bands  was  provided  by  the  Environmental 
Research  Institute  of  Michigan  (ERIM)  for  analysis  as  part  of  this  con- 
tract. A description  of  that  tape  is  paraphrased  below  from  a memo- 
randum provided  with  the  tape. 

A computer  compatible  data  tape  (CCT)  was  prepared  that  contains 
calibrated  images  from  Point  Mugu  and  Nellis  AFB.  Additional  calibrated 
images  will  be  prepared  in  the  future  but  they  have  not  been  converted 
from  high  density  digital  tape  (HDDT)  to  CCT  at  the  present  time.  The 
calibration  procedure  is  described  in  [1]  and  some  statistics  derived 
from  the  data  on  the  tape  are  included  in  [2]. 

The  nine  track  1600  BPI  data  tape  contains  9 files  (images).  The 
first  two  images  are  not  good.  Each  file  contains  a title  record  followed 
by  data  records.  Each  data  record  contains  one  scan  line  of  data.  Data 
are  put  onto  the  tape  with  four  9-bit  binary  data  values  per  (36-bit) 
word.  Data  values  are  limited  in  range  from  0 to  255  for  ease  of  pro- 
cessing on  a 32-bit  word  machine.  The  first  eight  (36-bit)  words  in  each 
data  record  are  header  information  for  the  scan  line.  The  first  (36-bit) 
word  is  a scan  line  counter.  Multispectral  data  are  channel  interleaved 
i.e.,  pixel  1 channel  1,  pixel  1 channel  2,  pixel  1 channel  3,  . . . 
pixel  646  channel  N. 

The  CCT  was  converted  to  32-bit  format  at  Auburn  for  ease  of  reading 
on  the  IBM370.  The  conversion  was  validated  by  comparing  selected  data 
values  from  the  converted  tape  with  values  provided  by  ERIM. 
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The  data  values  all  fall  in  the  range  0 to  255.  To  convert  these 

2 

digital  numbers  to  radiances  (yW/cm  • sr  • ym  for  the  2. 0-2. 6 ym 
channel  or  apparent  temperatures  (K  for  the  3. 0-4. 2,  3. 5-3. 9,  3. 9-4. 7, 

4. 5-5. 5,  5. 1-5. 7,  9.0-11.4  ym  channel),  they  have  to  be  scaled  as  follows: 
Calibrated  Data  = (digital  number )*MF  + AF 
The  images,  number  of  channels,  spectral  bands  in  order  as  they  are 
on  the  tapes,  number  of  scan  lines,  number  of  pixels  per  scan  line,  and 
MF  and  AF  for  each  channel  are  given  below. 

File  3,  NEVA!,  channels  2. 0-2. 6,  3. 0-4. 2,  4. 5-5. 5,  5. 1-5.7 
271  records  646  pixels,  2584  bytes  per  record 


MF  AF 


2.0  - 2.6 

3.0932 

-14.9566 

(approximate) 

3.0  - 4.2 

Bad  Data 

4.5  - 5.5 

.0316129 

283.956 

5.1  - 5.7 

.255415 

261.200 

(low  SNR  calibration 
very  uncertain) 

NEVA4,  channels  3.9-4. 7, 

3. 5-3. 9,  9.0-11 

.4 

215  records 

646  pixels. 

1938  bytes  per  record 

MF 

AF 

3.9  - 4.7 

.178320 

277.694 

3.5  - 3.9 

. 626448 

226.579 

9.0  - 11.4 

.0704379 

284.457 
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File  5,  NEVB,  channels  3. 9-4. 7,  3. 5-3. 9,  5.1-5. 7 


899  records  646  pixels. 

1938  bytes  per 

record 

MF 

AF 

3.9  - 4.7 

.085632 

277.476 

3.5  - 3.9 

.484893 

201.702 

5.1  - 5.7 

.517669 

199.105 

NEVC,  channels  3. 9-4. 7, 

3. 5-3. 9,  9.0-11 

.4 

898  records  646  pixels. 

1938  bytes  per 

record 

MF 

AF 

3.9  - 4.7 

.155668 

276.012 

3.5  - 3.9 

.504799 

224.004 

9.0  - 11.4 

.14172 

279.92 

File  8,  NEVD,  channels  3. 9-4. 7,  3. 5-3. 9,  9.0-11.4 

1541  records  646  pixels,  1938  bytes  per  record 


MF 


AF 


3.9  - 4.7 

.0966366 

267.812 

3.5  - 3.9 

.469730 

193.632 

9.0  - 11.4 

.127765 

269.867 

NEVE,  channels  3. 9-4. 7, 

3. 5-3. 9,  5. 1-5. 7 

179  record, 

, 646  pixels. 

1938  bytes  per  record 

MF 

AF 

3.9  - 4.7 

.127561 

268.445 

3.5  - 3.9 

. 644889 

173.319 

5.1  - 5.7 

.537182 

193.200 
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File  9,  NEVF,  channels  3. 9-4. 7,  3. 5-3. 9,  9.0-11.4 

899  records  646  pixels,  1938  bytes  per  record 
MF  AF 


3.9  - 

4.7 

.111736 

270.111 

3.5  - 

3.9 

.538573 

193.809 

9.0  - 

11.4 

.139151 

273.470 

All  of  the  calibrated  images  have  been  processed  to  some  degree. 

The  angular  resolution  of  the  scanner  is  2.5  mrad  and  the  data  are 
sampled  every  2.5  mrad  along  the  scan  line.  However,  the  scanner 
scans  60  lines/sec  and  the  aircraft  moves  202  ft/sec  so  that  there  is  some 
scan  line  overlap  in  data  collected  for  altitudes  above  1350  ft.  and  in 
the  data  collected  with  the  scanner  tilted  to  look  ahead  of  the  aircraft 
in  a 35°  depression  angle  configuration.  The  processing  that  has  been 
performed  on  the  original  data  is  one  of  line  averaging  yielding  an  along- 
track  ground  sample  distance  equal  to  the  along  track  ground  resolution. 
The  effective  ground  resolution  and  ground  sample  distance  for  each  image 
is  given  in  the  table  on  the  following  page. 

There  are  five  effects  in  the  imagery  and  in  the  data  that  should  be 
noted  in  any  further  processing  of  the  data.  (1)  NEVA!  and  NEVA4  are  not 
continuous  images.  From  an  entire  flight  of  several  miles  and  7000  scan 
lines,  about  45C  of  every  500  were  dropped  in  order  to  reduce  the  amount 
of  data  to  be  precessed.  Note  these  imagery  were  collected  to  determine 
the  magnitude  of  the  sunglint  and  such  sampling  is  adequate  with  the  air- 
craft making  a slow  level  turn  normal  to  the  suns  azimuth;  (2)  Some  of  the 
data  were  obtained  with  a linear  array  of  2 or  3 detectors.  Hence,  the 
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data  on  the  tape  are  not  in  spatial  registration.  This  should  be  corrected 
for  if  correlations  are  to  be  measured.  To  bring  the  data  into  registra- 
tion you  have  to  shift  the  data  as  follows: 

NEVA1  3. 0-4. 2 Shift  Left  2 Pixels 

4. 5- 5. 5 Shift  Right  2 Pixels 

i.e.,  registered  pixels  are 

(pi>u.)2>o_2.6’  ^pix^5.1-5.7’ 

(plx  *+2)3  Q_4  2,  (pix  x.-2)4  5_5  5 

All  Others  3. 9-4. 7 Shift  Right  2 Pixels 

3. 5- 3. 9 Shift  Left  2 Pixels 

i.e.,  registered  pixels  are 
(pixc) 

5. 1-5. 7 or  9.0-11.4 
(pix  ^"2 ) 3. 9_4. 7 

(pix  ^+2)3. 5-3.9 
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Data  Summary 


A/C  Alt. 
Above  Terrain 
(Ft) 

Scanner 
Depression 
(Deg. ) 

Cross-Track 

At  Nadir  (Ft. ) 

Along-Track 

At  Nadir  (Ft. ) 

Res 

GSD 

Res 

GSD 

NEVA1 

2000 

90 

5.0 

5.0 

5.0 

5.0 

NEVA4 

2000 

90 

5.0 

5.0 

5.0 

5.0 

NEVB 

1000 

35 

4.36 

4.36 

7.60 

7.60 

NEVC 

1000 

35 

4.36 

4.36 

7.60 

7.60 

NEVD 

1750 

90 

4.38 

4.38 

4.38 

4.38 

NEVE 

5000 

35 

12.5 

12.5 

38.00 

38.00 

NEVF 

1000 

35 

4.36 

4.36 

7.60 

7.60 

Date 

Time 

Heading 

Terrain  Type 

NEVA1 

3-7-78 

1130 

SW-SE 

Water 

NEVA4 

3-7-78 

1156 

SW-SE 

Water 

NEVB 

2-25-78 

1511 

E 

Mountains 

NEVC 

2-25-78 

1056 

U 

Desert 

NEVD 

2-25-78 

0914 

W 

Mountains 

NEVE 

2-25-78 

1424 

E 

Mountains 

NEVF 

2-25-78 

1034 

E 

Mountains 

(3)  Because  of  the  scanning  optics  the  image  of  the  detector  array  on  the 
ground  rotates  as  the  scan  mirror  rotates.  Directly  beneath  the  aircraft 
and  for  scan  angles  to  + 20  degrees  the  software  implemented  time  delays 
in  (2)  above  should  bring  the  data  into  spatial  registration 


to  within  a half  pixel.  In  order  to  achieve  registration  for 
the  larger  scan  angles  (to  ± 45°  degrees)  one  would  have  to 
shift  pixels  from  different  lines,  but  this  is  probably  not 
a worthwhile  exercise  because  the  along  track  ground  resolu- 
tion begins  to  exceed  the  along-track  ground  sample  distance 

2 

along  the  scan  line  by  a factor  of  sec  9.  (4)  In  the  tilted 

scanner  data  there  is  some  vignetting  by  the  aircraft  skin 
at  the  edges  of  the  scans.  Finally,  (5)  NEVC  is  an  image 
of  desert  terrain  at  Nellis  and  the  image  is  actually  composed 
of  two  discontinuous  segments  of  a larger  image. 


[1]  Beard,  J.,  Maxwell,  J.  R.  and  Spellicy,  R. , Statistical 
Analysis  of  Terrain  Background  Measurements  Data,  ERIM 
Report  No.  120500-12-F,  Environmental  Research  Institute 
of  Michigan,  Ann  Arbor,  March  1977. 

[2]  Maxwell,  J.  Robert,  "Statistical  Analyses  of  Selected 
Terrain  and  Water  Background  Measurements  Data," 

ERIM,  Report  No.  132300-1-F,  July  1978. 
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SECTION  IV.  SOFTWARE  DESCRIPTION  FOR  DATA  ANALYSIS 

Development  of  software  for  this  contract  has  resulted  in  the  writing 
of  three  main  programs.  These  are  Plotit.  Area  and  Shade.  Plotit  calcu- 
lates amplitude  statistics  and  generates  amplitude  histograms  in  terms  of 
2 

radiance  (W/M-SR).  Area  defines  blocks  of  connected  “picture"  elements 
whose  radiance  level  is  above  a threshold.  It  also  generates  a histogram 
of  number  of  elements  in  a block  vs.  their  relative  occurrence.  Shade 
generates  a "picture"  from  the  data  sets  by  linking  areas  corresponding 
to  data  values  greater  than  a threshold. 

To  understand  Plotit  we  must  first  look  at  the  subroutines  it  uses. 
These  are  Rad,  Quint,  and  Lion.  Rad  first  converts  an  integer  data  value 
between  0 and  255  to  an  equivalent  blackbody  temperature  using  the  equation 
T = MF  (Data  Value)  + AF  where  MF  and  AF  are  given.  Then  using  the 
equivalent  temperature.  Rad  numerically  integrates  Planck's  equation  over 
the  spectral  band  to  determine  the  radiance  level  equivalent  to  the  input 
integer  data  value.  Subroutine  Quint  does  the  inverse  of  It  finds 

the  nearest  integer  equivalent  to  a radiance  level  wnich  is  inputed  to  it. 
Subroutine  Lion  generates  axis  and  labeling  for  the  histogram.  It  then 
scales  and  plots  the  data  presented  to  it  in  two  1-D  arrays,  x and  Y. 

The  main  program  of  Plotit  first  calculates  values  for  the  two  arrays  to 
be  used  for  the  amplitude  histogram.  X(I)  is  equal  to  the  radiance 
corresponding  to  the  integer  1-1.  Y ( I ) is  equal  to  the  number  of  data 
elements  equal  to  1-1  divided  by  the  total  number  of  data  elements.  For 
X(I)  and  Y(I)  I ranges  from  1 to  256.  The  main  program  then  calls  Lion. 
From  the  arrays  X ( I ) and  Y ( I ) the  main  program  next  calculates  unbiased 
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and  consistent  estimates  for  the  mean  (u),  variance  (a2),  and  standard 
deviation  (a)  by  use  of  the  equations  below. 


256 

X(I)  Y(I) 

(9) 

1=1 

256 

'2  N 

a = fPT 

[(][]  cx(i)]2  y ( i ) ) . ;2] 

(10) 

1=1 

a = 

(ID 

where  N equals  the  total  number  of  pixels  in  a given  image.  Finally 
P^otit  prints  out  the  integer  quantization  values  equivalent  to  radiance 
levels  of  u + .25na  for  n equal  0 to  16.  These  values  are  used  as  thresholds 
in  programs  Area  and  Shade. 

The  programs  Shade  and  Area  consider  the  "scanned  picture"  nature 
of  the  data  sets.  They  consider  each  data  element  as  a square  in  a grid. 


Scan 

line 

or  2 

I 

Record  • 

I 

I 


M 


For  example  the  second  element  of  scan  line  M is  the  square  labled  A in 
the  figure  above. 


Element 
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The  program  Shade  generates  a 2-0  plot  by  darkening  squares  corre- 
sponding to  data  points  above  a selectable  threshold  and  leaving  white 


i 


all  other  blocks.  The  Program  also  draws  a rectangle  delimiting  the 
bounds  of  the  grid.  To  save  CPU  time  the  program  darkens  continuous 
blocks  along  a scan  line  at  the  same  time. 

Program  Area  is  the  most  complicated  of  the  three  programs.  Its 
operation  is  best  explained  in  terms  of  an  example. 


m 


In  the  figure  above,  the  darkened  blocks  represent  pixels  above  a given 
threshold.  Area  would  determine  that  the  data  set  this  figure  represents 
has  two  blocks  of  one,  two  blocks  of  three,  one  block  of  four  and  three 
blocks  of  eight.  Thus,  there  are  eight  blocks  in  this  example  data  set. 
Notice  that  the  squares  encircled  are  one  block  of  four.  Area  calculates 
the  value  of  two  arrays  X(I)  and  Y(I).  X(I)  = I for  I equal  1 to  1000. 
Y(I)  = number  of  blocks  of  I divided  by  the  total  number  of  blocks  for  I 
equal  1 to  999.  Y(1000)  = the  number  of  blocks  with  1000  or  more  elements 

divided  by  the  total  number  of  blocks.  Using  these  arrays,  a histogram 
is  generated  with  a subroutine  similar  to  Lion  in  Plotit.  It  is  plotted 
for  I equal  one  to  200.  Nonzero  values  of  Y(I)  for  201  <_  I <_  1000  are 
printed  out  for  evaluation.  The  option  to  print  all  histogram  values 
is  available  also. 
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The  above  criteria  for  plotting  only  the  first  200  area  histogram 
values  and  for  lumping  all  areas  above  1000  were  based  on  results  obtained 
from  analyzing  the  data.  In  fact  the  area  histograms  are  heavily  biased 
toward  areas  of  twenty  or  fewer  pixels.  The  point  is  that  these  plot 
parameters  can  be  adjusted  to  best  show  scene  dependent  area  statistics. 

The  computer  generated  plots  were  produced  on  a Versatec  electro- 
static plotter,  which  when  properly  adjusted  should  have  the  capability 
of  generating  gray  scale  plots.  Software  for  gray  scale  plots  was 
generated  and  tested  using  ten  shades  of  gray  as  one  inch  bars  across  the 
plotter  width.  Difficulty  was  encountered  with  the  darkest  gray  levels. 
Additionally  the  gray  scale  plot  was  inefficient  because  of  the  required 
mechanism  for  addressing  the  Versatec. 

This  program  was  abandoned  because  of  the  above  mentioned  problems 
and  because  other  methods  for  gray  scale  image  production  already  exist. 
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SECTION  V.  RESULTS  OF  ANALYSIS 


Results  of  the  analysis  consist  of  the  following 

5.1  Data  Summary  Printout  - For  each  spectral  band  and  data  file, 
the  following  information  is  printed:  MEAN,  STANDARD  DEVIATION, 
VARIANCE,  NUMBER  OF  PIXELS,  MINIMUM  RADIANCE  (Corresponding  to 
bin  number  0),  MAXIMUM  RADIANCE  (Corresponding  to  bin  number 
255),  and  BIN  NUMBERS  for  selected  fractional  standard  devia- 
tions above  the  mean. 

5.2  Area  Histogram  Printout  - For  each  spectral  band,  file  and 
threshold,  the  areas  above  200  pixels  in  size  are  printed  out 
along  with  a total  count  of  discrete  areas  for  the  scene  and 
spectral  band. 

5.3  Amplitude  Histogram  Plots  - For  each  file  and  spectral  band 
an  amplitude  histogram  based  on  all  pixels  is  computed  and 
plotted.  Conversion  of  the  data  to  radiances  is  accomplished 
prior  to  histogram  computation. 

5.4  Shade  Plots  - For  each  file,  spectral  band  and  threshold,  a 
two-dimensional  plot  is  generated  showing  only  those  pixels 
above  the  threshold.  The  threshold  is  chosen  as  some  fractional 
standard  deviation  above  the  mean.  Example  plots  are  included 
in  this  report.  Because  of  the  large  shaded  areas,  these  plots 
had  to  be  reproduced  photographically. 

5.5  Area  Histogram  Plots  - For  each  file,  spectral  band  and  thresh- 
old as  generated  in  5.4,  a histogram  of  areas  is  computed  and 
plotted.  Examples  included  verify  the  bias  toward  small  areas. 


! 
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The  following  pages  illustrate  results  of  the  various  computations. 
Identifying  information  is  given  for  each  figure. 
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INFORMATION  FOR  INTERPRETING  LABELS 


FILE 

SPECTRAL  BAND  NUMBER 

WAVELENGTH  (MICRONS) 

NAME 

5 

1 

3.9 

- 4.7 

NEVB 

5 

2 

3.5 

- 3.9 

NEVB 

5 

3 

5.1 

- 5.7 

NEVB 

6 

1 

3.9 

- 4.7 

NEVE 

6 

2 

3.5 

- 3.9 

NEVE 

6 

3 

5.1 

- 5.7 

NEVE 

7 

1 

3.9 

- 4.7 

NEVC 

7 

2 

3.5 

- 3.9 

NEVC 

7 

3 

9.0 

-11.4 

NEVC 

8 

1 

3.9 

- 4.7 

NEVD 

8 

2 

3.5 

- 3.9 

NEVD 

8 

3 

9.0 

-11.4 

NEVD 

9 

1 

3.9 

- 4.7 

NEVF 

9 

2 

3.5 

- 3.9 

NEVF 

9 

3 

9.0 

-11.4 

NEVF 

DATA  SUMMARY 


FILE  5 


SPECTRAL  BAND  1 

Mean  = 5.76  x 10'  1 W M'2  SR"1 
Standard  Deviation  = 9.96  x 10"2  W M~2  SR-1 
Variance  = 9.92  x 10' 5 W2  M"4  SR'2 
Number  of  Pixels  = 580108 

- 1 -2-1 

Maximum  Radiance  = 5.96  x 10  W M SR 

-1  -?  -1 

Minimum  Radiance  = 5.53  x 10  W M SR 

SPECTRAL  BAND  2 

Mean  = 5.28  x 10" 2 W M"2  SR"1 

Standard  Deviation  = 1.36  x 10" 3 W M'2  SR'1 

Variance  = 1.85  x 10' 6 W2  M"4  SR'2 

Number  of  Pixels  = 580108 

Maximum  Radiance  = 5.80  x 10'2  W M~2  SR'1 

Minimum  Radiance  = 3.60  x 10’ 2 W M'2  SR’1 

SPECTRAL  BAND  3 

Mean  =1.67  W M'2  SR'1 

Standard  Deveiation  = 2.24  x 10' 2 W M'2  SR 

Variance  = 5.01  x 10'^  W2  M’4  SR’2 

Number  of  Pixels  = 580108 

Maximum  Radiance  =1.95  W M’2  SR’1 

Minimum  Radiance  =1.17  W M’2  SR'1 
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DATA  SUMMARY  FILE  6 
SPECTRAL  BAND  1 

Mean  = 5.76  x 10' 1 W M"2  SR"1 
Standard  Deviation  = 1.15  x 10' 2 W M~2  SR'1 
Variance  = 1.32  x 10' 4 W2  M'4  SR"2 
Number  of  Pixels  = 114988 


Maximum  Radiance  = 
Minimum  Radiance  = 

Mean  = 5.43  x 10'  2 
Standard  Deviation 
Variance  = 2.50  x 
Number  of  Pixels  = 
Maximum  Radiance  = 
Minimum  Radiance  = 

Mean  =1.64 
Standard  Deviation 
Variance  = 2.04  x 

[ I 

Number  of  Pixels  = 
Maximum  Radiance  = 
Minimum  Radiance  = 


6.00  x 10"1  W M"2  SR'1 

5.35  x 10'1  W M"2  SR'1 

SPECTRAL  BAND  2 

-2  -1 
W M c SR  1 

= 1.58  x 10" 3 W M'2  SR'1 
0"6  W2  M"4  SR"2 
114988 

6.07  x 10' 2 W M"2  SR"1 

3.09  x 10"2  W M'2  SR'1 


SPECTRAL  BAND  3 


W M" 

2 SR"1 

= 1.43  x 

10"2 

W 1 

m"2 

SR 

vf 

1 

O 

W2  M'4 

SR'2 

114488 

1.95 

W M"2 

SR' 

■1 

1.14 

W M"2 

SR' 

■1 
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DATA  SUMMARY  FILE  7 
SPECTRAL  BAND  1 


Mean  = 5.91  x 10'1 

W M 

-2  -1 

1 ^ SR 

Standard  Deviation 

= 6.26 

x 10"  3 

W M 

r2 

SR 

Variance  = 3.92  x 

10"5 

W2  M"4 

SR"2 

Number  of  Pixels  = 

579362 

Maximum  Radiance  = 

6.29  x 

10"  1 

W M"2 

SR* 

■1 

Minimum  Radiance  = 

5.50  x 

10"1 

W M'2 

SR* 

■1 

SPECTRAL  BAND  2 

Mean  = 5.51  x 10" 2 W M"2  SR"1 

Standard  Deviation  = 1 . 22  x 10"-^  W M"2  SR 

Variance  = 1.48  x 10" 6 W2  M"4  SR'2 

Number  of  Pixels  = 579462 

Maximum  Radiance  = 6.29  x 10’ 2 W M"2  SR"1 

Minimum  Radiance  = 3.99  x 10" 2 W M"2  SR'l 

SPECTRAL  BAND  3 

Mean  = 5.71  x 10'10  W M"2  SR'1 

Standard  Deviation  = 6.41  x 10’  2 W M 2 SR 

Variance  = 4.10  x TO"'1  W2  M’4  SR"2 

Number  of  Pixels  = 579462 

Maximum  Radiance  =6.08  W M”2  SR"1 

Minimum  Radiance  =5.38  W M"2  SR"1 
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DATA  SUMMARY  FILE  8 


SPECTRAL  BAND  1 

Mean  = 5.66  x 10"1  W M"2  SR'1 
Standard  Deviation  = 8.93  x 10' 3 W M’2  SR"1 
Variance  = 7.98  x 10"5  W2  M"4  SR'2 
Number  of  Pixels  = 994840 
Maximum  Radiance  = 5.83  x 10" ^ 

Minimum  Radiance  = 5.34  x lO"1- 

SPECTRAL  BAND  2 

Mean  = 5.21  x 10' 2 W M"2  SR"1 

Standard  Deviation  = 1.72  x 10"  3 W M"2  SR"1 

Variance  = 2.47  x 10'6  W2  M"4  SR'2 

Number  of  Pixels  = 994840 

Maximum  Radiance  = 5.59  x 10" 2 W M"2  SR"1 

Minimum  Radiance  = 3.45  x 10"2  W M"2  SR'1 

SPECTRAL  BAND  3 

Mean  =5.48  W M"2  SR"1 

Standard  Deviation  = 1.06  x 10"1  W M"2  SR"1 
Variance  = 1.13  x 10"2  W2  M*4  SR"2 
Number  of  Pixels  = 994840 

? - 1 

Maximum  Radiance  = 5.81  W M SR 

Minimum  Radiance  =5.19  W M"2  SR"1 


DATA  SUMMARY  FILE  9 


SPECTRAL  BAND  1 

Mean  = 5.70  x 10" 1 W M'2  SR'1 

Standard  Deviation  - 1.02  x 10' 2 W M'2  SR'1 

Variance  = 1.04  x 10'  4 W2  M"4  SR'2 

Number  of  Pixels  = 580108 

Maximum  Radiance  = 5.95  x 10~1  W M~^  SR~ 

Minimum  Radiance  = 5.38  x 10’^-  W M"2  SR"1 

SPECTRAL  BAND  2 

Mean  = 5.27  x 10'2  W M'2  SR'1 
Standard  Deviation  * 1.63  x 10*^  W M"2  SR'1 

Variance  = 2.67  x 10'  6 W2  M'4  SR'2 

Number  of  Pixels  = 580108 

_ o _2-l 

Maximum  Radiance  = 5.90  x 10  W M SR 

Minimum  Radiance  = 3.45  x 10" 2 W M"2  SR"1 

SPECTRAL  BAND  3 

Mean  =5.52  W M'2  SR'1 

Standard  Deviation  = 1.19  x 10"1  W M"2  SR"1 

Variance  = 1.41  x 10'2  W2  M'4  SR'2 

Number  of  Pixels  = 580108 

Maximum  Radiance  =5.95  W M~2  SR"1 

-?  -1 

Minimum  Radiance  =5.26  W M SR 
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SUMMARY  OF  DATA  NOT  PLOTTED  ON  AREA  HISTOGRAMS 


OF  FILE  5,  SPECTRAL  BAND  1 
THRESHOLD  = y + a 


NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10 


-3 


205 

213 

222 

256 

495 

529 

805 

977 

1000  or  Greater 


1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

1.6 

19.0 


TOTAL  NUMBER  OF  BLOCKS  = 644 


THRESHOLD  = y + 1.5a 


NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10 


-3 


206 

210 

238 

253 

313 

393 

435 

441 

565 

661 

1000  or  Greater 


1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

15.0 


TOTAL  NUMBER  OF  BLOCKS  = 786 


THRESHOLD  = y + 2a 
TOTAL  NUMBER  OF  BLOCKS  = 325 
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SUMMARY  OF  DATA  NOT  PLOTTED  ON  AREA  HISTOGRAMS 
OF  FILE  5,  SPECTRAL  BAND  2 
THRESHOLD  = u + a 

NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10‘4 


213  2.5 

281  2.5 

336  2.5 

354  2.5 

397  2.5 

486  2.5 

565  2.5 

874  2.5 

1000  or  Greater  32.0 


TOTAL  NUMBER  OF  BLOCKS  = 4070 

THRESHOLD  = y + 1.5a 

NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10'4 


206  2.9 

242  2.9 

334  2.9 

872  2.9 

1000  or  Greater  18.0 


TOTAL  NUMBER  OF  BLOCKS  = 3404 

THRESHOLD  = v + 2a 
TOTAL  NUMBER  OF  BLOCKS  = 426 
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SUMMARY  OF  DATA  NOT  PLOTTED  ON  AREA  HISTOGRAMS 


OF  FILE  5,  SPECTRAL  BAND  3 
THRESHOLD  = y + 1.5a 

NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10-5 
288  7.2 

TOTAL  NUMBER  OF  BLOCKS  = 13796 

THRESHOLD  = y + 2 a 
TOTAL  NUMBER  OF  BLOCKS  = 6886 

• THRESHOLD  = y + 2.5a 
TOTAL  NUMBER  OF  BLOCKS  = 1441 

THRESHOLD  = y + 3a 
TOTAL  NUMBER  OF  BLOCKS  = 161 
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SUMMARY  OF  DATA  NOT  PLOTTED  ON  AREA  HISTOGRAM 


OF  FILE  7,  SPECTRAL  BAND  1 
THRESHOLD  = u + o 


NUMBER  OF  PIXELS  IN  A BLOCK 

209 

216 

213 

217 

231 

235 

244 

253 

266 

284 

297 

306 

326 

347 

354 

364 

438 

443 

454 

477 

590 

604 

712 

904 

934 

961 

962 

1000  or  Greater 


RELATIVE  OCCURRENCE  X 10'4 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

3.5 

28.0 


TOTAL  NUMBER  OF  BLOCKS  = 2837 

THRESHOLD  = p + 1 . 5c 


NUMBER  OF  PIXELS  IN  A BLOCK 

225 

437 

567 

863 

TOTAL  NUMBER  OF  BLOCKS  = 1030 


RELATIVE  OCCURRENCE  X 10"4 

9.7 

9.7 

9.7 

9.7 


THRESHOLD  = u + 2c 


TOTAL  NUMBER  OF  BLOCKS  = 14 
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SUMMARY  OF  DATA  NOT  PLOTTED  ON  AREA  HISTOGRAMS 
OF  FILE  7,  SPECTRAL  BAND  2 
THRESHOLD  = y + a 

NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10‘4 


203  2.5 

217  2.5 

255  2.5 

564  2.5 

666  2.5 

933  2.5 

1000  or  Greater  20.0 


TOTAL  NUMBER  OF  BLOCKS  = 3950 

THRESHOLD  = y + 1.5a 

NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10 


206  5.3 

334  5.3 

356  5.3 

466  5.3 

543  5.3 

560  5.3 

750  5.3 

1000  or  Greater  16.0 


TOTAL  NUMBER  OF  BLOCKS  = 1893 

THRESHOLD  = y + 2a 
TOTAL  NUMBER  OF  BLOCKS  = 376 

THRESHOLD  = y + 2.5a 
TOTAL  NUMBER  OF  BLOCKS  = 2 


i 
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SUMMARY  OF  DATA  NOT  PLOTTED  ON  AREA  HISTOGRAMS 
OF  FILE  7,  SPECTRAL  BAND  3 
THRESHOLD  = y + a 

NUMBER  OF  PIXELS  IN  A BLOCK  RELATIVE  OCCURRENCE  X 10'4 


226  2.6 

245  2.6 

253  2.6 

254  2.6 

301  2.6 

308  2.6 

397  2.6 

664  2.6 

665  2.6 

731  2.6 

741  2.6 

1000  or  Greater  7.8 


TOTAL  NUMBER  OF  BLOCKS  = 3858 

THRESHOLD  = y + 1.5a 
TOTAL  NUMBER  OF  BLOCKS  = 337 

THRESHOLD  = y + 2a 
TOTAL  NUMBER  OF  BLOCKS  = 4 
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AMPLITUDE  HISTOGRAM  FILE  5 SPECTRAL  BAND  2 


RRO I flNCE 


AMPLITUDE  HISTOGRAM  FILE  5 SPECTRAL  BAND  3 


i 


oi  ,ro 


80  J0  90  J0 

33N3yyn33o 


TlO  • 0 

3AI1U1 


a 3 


00  -0° 
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.08  0.10  0.12  0.14  0.16  0.18  0.20  0.22 

RRDIRNCE  W M"2  SR"1  *10 


AMPLITUDE  HISTOGRAM  FILE  6 SPECTRAL  BAND  1 


RRDIflNCE  W M-2  SR" 


AMPLITUDE  HISTOGRAM  FILE  7 SPECTRAL  BAND  2 


RADIANCE 


AMPLITUDE  HISTOGRAM  FILE  7 SPECTRAL  BAND  3 


ROD  I R N C E 


AMPLITUDE  HISTOGRAM  FILE  8 SPECTRAL  BAND  1 


RADIANCE  W M'2  SR- 


AMPLITUDE  HISTOGRAM  FILE  8 SPECTRAL  BAND  3 


RADIANCE  N M~2  SR” 


AMPLITUDE  HISTOGRAM  FILE  9 SPECTRAL  BAND  1 


5.50  5.60  5.70 

RADIANCE  N M'2  SR"1 


SHADE  PLOT  FILE  5 SPECTRAL  BAND  1 
THRESHOLD  = p + 2a 

W KVB  — t lOUilTAINS  (3 .9-4. 7 UICROHO) 


SHADE  PLOT  FILE  5 SPECTRAL  BAND  2 
THRESHOLD  = p + 1.55 
NEVB-MOUNTAINS  (3. 5-3. 9 MICRONS) 


SHADE  PLOT  FILE  5 SPECTRAL  BAND  3 
THRESHOLD  =0+1.5 5 

HEVB -MOUNTAINS  (5. 1-5. 7 MICRONS) 


SHADE  PLOT  FILE  5 SPECTRAL  BAND  3 
THRESHOi  D = fi  + 2.0a 

MF.VB-MOUNTAINf!  (5.1-5. 7 MTCRONS) 


SHADE  PLOT  FILE  5 SPECTRAL  BAND  3 
THRESHOLD  = ji  + 2.53 

HEVB-ilOUNTAIHC  (5. 1-5. 7 1IICRONS) 


SHADE  PLOT  FILE  5,  SPECTRAL  BAND  3 
THRESHOLD  = v + 3.0o 

NEVB -MOUNTAINS  (5.1-5. 7 MICRONS) 


[CR< 
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SHADE  PLOT  FILE  7 SPECTRAL  BAND  1 
THRESHOLD  = + 1.5a 

W I£VC -DESERT  (3. 9-4. 7 1IICRONS) 


SHADE  PLOT  FILE  7 SPECTRAL  BAUD  1 
THRESHOLD  = g + 2a 
NEVC-DESEUT  (3.9-4. 7 MICRONS) 


SHADE  PLOT  FILE  7 SPECTRAL  BAND  2 
THRESHOLD  = » + a 

N LVC- DESERT  (3. 5-3. 9 MICRONS) 


SHADE  PLOT  FILE  7 SPECTRAL  BAND  2 
THRESHOLD  = p + 2 a 
NEVC- DESERT  (3. 5-3. 9 MICRONS) 


SHADE  PLOT  FILE  7 SPECTRAL  BAND  2 
THRESHOLD  = y + 2.55 
UEVC-DESERT  (3.5-3. 9 MICRONS) 


PLOT  FILE  7 SPECTRAL  BAND 
THRESHOLD  =0+5 

DESERT  (9.0-11.4  MICRONS) 


SHADE  PLOT  FILE  7 SPECTRAL  BAUD  3 
THRESHOLD  = \>  + 1 . 5o 
WCVC-DESERT  (9.0-11.4  HICROtJS) 


AREA  HISTOGRAM 


AREA  HISTOGRAM 


.00  40.00  80.00  120.00  160.00  200.00  240.00  280.00 

NUMB E H OF  PIXELS  IN  BLOCK 


AREA  HISTOGRAM 


AREA  HISTOGRAM 


AREA  HISTOGRAM 


.00  110.00  80.00  120.00  160.00  200.00  240.00  280.00 


AREA  HISTOGRAM 


AREA  HISTOGRAM 


80.00  120.00  160.00  200.00  240.00  280.00 

NUMBER  OF  PIXELS  IN  BLOCK 


AREA  HISTOGRAM 


80.00  120.00  160.00  200.00  240.00  280.00 

NUMBER  OF  PIXELS  IN  BLOCK 


SECTION  VI.  SUMMARY  AND  CONCLUSIONS 


A brief  theoretical  development  was  presented  to  establish  a basis 
for  using  background  statistics  in  a discrimination  model.  The  develop- 
ment was  limited  to  threshold  detection  and  did  not  include  detailed 
spectral,  spatial  or  temporal  dependence.  An  advanced  model  for  dis- 
crimination would  probably  include  these  features  as  well. 

The  major  emphasis  for  the  contract  was  the  development  of  statistics 
from  empirical,  background,  banded,  radiometric  data  which  shows  thresh- 
old and  area  features.  The  analysis  was  primarily  for  individual  scenes 
and  was  not  to  duplicate  existing  methods  for  inter-spectral  band  correla- 
tion. Thus  the  results  serve  as  one  specific  class  of  analysis  methods 
which  should  be  used  together  with  other  methods  in  a discrimination 
model . 

Application  of  the  analysis  to  more  suitable  background  scenarios 
should  provide  useful  inputs  for  infrared  warning  receiver  false  alarm 
analysis,  when  used  in  conjunction  with  realistic  target  signature 
measurements. 
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APPENDIX 
Software  Listing 


PLOTIT  GENERATES  AMPLITUDE  HISTOGRAMS  AND  CALCULATES  AMPLITUDE 
STATISTICS  FROM  DIGITIZED  MULTI  SPECIAL  SCANNEP  DATA  STORED  SCAN 
LINE  PER  RECORD  ON  MAGNETIC  TAPE-  PLOTIT  IS  SET  U»  TO  READ  TAPES 
ON  WHICH  DATA  VALUES  FOR  THREE  SPECTRAL  “ENDS  HAVE  BEEN 
INTERLEAVED.  THIS  PROGRAM  ASSUMES  THAT  THERE  ARE  1938  DATA 
VALUES  PER  RECORO  and  That  VALUES  RELONG  TC  EACH  SPECTRAL 

SAND.  IT  ALSO  ASSUUES  T*AT  THE  MAGNITUDE  CF  DATA  VALUES  RANGES 
FROM  0 TO  255.  JOB  CCNTPCL  LANGUAGE  MUST  3E  INCLUDED  TO  IDENTIFY 
UNIT  -2  AS  T A°  E CONTAINING  FILE  TO  BE  USED. 

DIMENSION  XI258) ,Y(255! . I D ( 19351 
REAL  MF 

INTEGER  SB.QU.OUSD 

,***««**  PARAMETERS  WHICH  CHANGE  WITH  DATA  SETS  ***************** 
NSR=898 
F I LE=5 
SB  s I 

AF*277.A7 6 
MF*. 085632 
UM I N=  3 . 5 
UMAX  = 3.  9 

**************************=**■,*************  *****************  r.  -•****. 
NBR  = NUMSER  OF  0RCCr'’'S  IN  FILE 

FILE=FIL£  F°Om  w-j ICH  Data  IS  tskFN  USED  ONLY  FOR  LABELING 
S8=SPECTRAL  RAND  RR  interest  1.2.0°  3 

AF  AND  MF  A°E  CAL  A B-R  AT  I C\'  CONSTANTS  SU?»LYED  WITH  TAPES 
UMIN=LOWER  LIMIT  3C  SPECTRAL  BAND  (MICRONS) 

UMAX=U°PE°  LIMIT  CF  SPECTRAL  RAND  (“ICRCNSI 


NP  S=  256 
VALUE-SB 
N=0 
U=0. 

X2  = 0. 

C ****** ******  S°T  U°  Y ARRAV  a***************************************** 

DO  5 1=1,256 
5 Y(II=0. 

DO  10  K= 1 , NBR 
RE  ADI  2 I ID 
DO  10  I =S  B , 1938, 3 
N = N*1 
J=IO( I 1+1 
10  Y(J)=Y(J)+1 

DO  20  K = l , 2 56 
IFIYIKI.EO.O.)  GO  TO  20 
Y(K)=Y(KI/N 

20  CONTINUE 

C************  SET  UP  X ARRAY  ************************************** **** 

DO  21  1=1,256 
A*  I - 1 

CALL  RADIA.AF.MF.UMIN.UMAX.R) 

21  X(I)=R 

c************  PLOT  HISTAGPAM  ****************************************** 

CALL  LION! X.Y.NPS, VALUE, FILE) 

C************  CALCULATE  MEAN,  VARIANCE,  AND  STANDARD  DEVIATION********* 
DO  22  1=1,256 
U=X(  I )*Y( I) *U 

22  X2=YII)*X( I )**2*X2 
S=( N/ ( N— 1 I ) *( X2-U**2 ) 


THIS  PA G£  IS  BESI  QUALITY  PKACIlCA&ui 
FROM  OOPY  PURiilSHED  10  DOC 


23 


2f> 
30 
35 
40 
u s 

50 


S0=S0RT(S) 

F I NO  INTEGEB^OUNTIZATION  LEVEL  CORRESPONDING  TC  THE  mean 

CALL  OUNTIU.X^OU)1 ' 

WRITEI6.30)  OU 

PINO  OUNTIZ4TION  LEVEL  CHRP  ESDONOTNG  TO  ME AN* ( . 25  I ( A H SO  I 
RAOIANCE  LEVEL  FOR  A=1  TO  40  **  = »»*»*•*»*»*•*«**»•*******<■******** 
DC  23  K=l,40 
A=K«.25 
R I N=U*SD*  A 

CALL  OUNTI PIN, X. OUSDJ 
WRITE16.40)  A , OUSD 
CONTINUE 

WRITE<6,50)  FILE 
WRITE(6,45I  UMAX,UMIN,MF,AF 
WRITF(6,26)  N 
WRITE(6,35IU,S.SD 

FORMAT!  'O'  , 'NUMBER  OF  POINTS  = ',16  ) 

FORMAT! '0* ,' ME AN=  ',131 

FORMAT! 'O' , 'MEAN  = '.Ell. 4,'  VARIANCE  = 'rE11.4,'  SO  = •, Ell-4) 
FORMAT! *0’ , 'MEAN  +',F4.2.'SD  = ’,13) 

FORMAT! 'O'  ,'UMATr  '.F4.1,'  UMIN  = ',F4.1.'  MP=  *,F9.7,'  A F=  ',F9.31 
FORMAT! '0* , '*»**»***  FILE  = *tF4.1»'  =*«»**=*•) 

STOP 

END 


C 
c 
c 
c 

SUBROUTINE  0UNT(rin,X,0! 

C******«******QUINI  FINOi  The  3UNT 1 2 A T I ON  level. 
C INOUT  PAOIANCE  LEVEL.  RIN.  IT  USES  ' 

C IN  MAIN  PROGRAM. 

DIMENSION  X ( 25  3 » 

INTEGER  0 
S»l. 

0=300 

00  5 1=1,256 

01  F=AB$IR IN-X! T ) I 
IFIDIF.GT.Sl  GO  TO  5 
S = 0IF 
0=1-1 

5 CONTINUE 
RETURN 
ENO 
C 
C 
C 

c 

SUBROUTINE  RAOIA.AF.MF.UMIN.UMAX.P) 
C*****«»***«»*S AO  CALCULATES  THE  RADIANCE  LEVEL. 

C THE  INTEGER  DATA  VALUE,  A. 

REAL  MF 
T=A*MF+AF 
P 1=3 .1415927 
Cl=3.7413E-08/PI 
C2=l .4388 
R = 0. 


= ,CC3R PS°ON0I NG  TO  AN 
HE  XARR.AY  CALCULATED 


,CCRR ESRONOING  TC 


rUCUCABLM 


TO  DDQ 
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U=UMIN+.005 

W*  Cl/I <EXP(C2/(U*T) >-l.)*U**5) 

»=R+w». 01 

U»U+.0l 

IF ( U.LT.UMAXI  GO  TO  5 

RETURN 

END 


SUBROUTINE  LI  ON (X.Y.NPS. VALUE. FILE) 

C*********«*«»L ION  GENERATES  ALL  A X I S , SC AL I NO , LABELING  AND  PLOT  I NG 


C NECESSARY  FOR  MAKING  HISTOGRAMS  •>.!»•.•»«»«*>»»»»«>•* 

0IMENS1CN  X ( 2 5 A ) . YI258) 

CALL  PLOYSIO.O.O) 

CALL  PLOTID.,1.,2) 

CALL  NUUQS°(C.,R.,.35.FIL=.0.,-l) 

CALL  NIJM3ER  ( 0.0.9.  5.  - 35,  VALUE.O.  0,-11 
CALL  PLOT ( 3. . I .25 ,-3 ) 

CALL  SCALE ( Y, 5. ,NPS. I ) 

- CALL  SCALE ( X, 7.0.NPS, 1 ) 

CALL  AXIS(0.,0., 'RELATIVE  OCCURRENCE' , 1 9 , R . , °C.»Y(NPSwl ) , 

9 YINPS+2)  ) > 

C»****«»t*****XAX I S IS  AMODIFIED  VERSION  OF  THE  VEGSATEC  SUBROUTINE 
C AXIS.  IT  ALLOWS  SUPERSCRIPTS  TC  A^EAR  IN  AXIS  LABEL.  »*** 

CALL  XAX I S ( 0. , 0. , ' RAD  I ANCc  W M -2  SR  -1  ' . -2 A , 7 . , 0 . , X I NP S ♦ 1) , 
aX(N°S+2) ) 

CALL  LINE! X.Y.NPS, 1.0,0 ) 

CALL  PLOT! 10. ,-1.25,-31 
CALL  PLOTIO.,1.,2) 

CALL  PL0T(C.,0.,*9O9l 

RETURN 

END 


OUTPUT  JOB  8895  9V  GOGGANS 


167  CARDS  READ  3 PAGES  PRINTED 


0 CA°DS  PJNCMED 


THIS  PA 

F u nacn. 
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I 


c shade  generates  a 2- o riot  from  digitized  multispectral  scanner 

C DATA  STCPED  CN  MAGNETIC  tape.  SuAOE  DARKENS  SQUARES  CCRR ESPQNOING 

C TO  DATA  VALUES  OVER  A SELECTABLE  THRESHOLD  AND  LEAVES  WHITE  ALL 

C OTHER  ARES.  SHADE  MAKES  THE  SA-E  ASSUMPTIONS  ABOUT  THE  DATA  AS 

C PROGRAM  PLCTIT.  JC  3 CCNT  = CL  LANGUAGE  MUST  3E  INCLUOEO  TO  IDENTIFY 

C UNIT  2 AS  TAPE  CONTAINING  FILE  TC  3 E USED.  A °LOT  PARAMETER  CARO 

C MUST  AL SC  ? E INCLUDED  TC  SET  °LCT ER  UNHS  = CNE  UNIT  PER  .01  INCH. 

INTEGER  THRES 

integer  oat  a< 667  i , s.c  ,a,t,sb,pre 

0 I MENS ICN  10(1933) 

DIMENSION  XI4I  ,Y( 4|  ,IPAT9t 16) 

DATA  IRAT9/16*ZFECF/ 

PARAMETERS  WHICH  CHANGE  WITH  OATA  SETS  ***************** 

M = 893 
THRES=169 
SB  = 3 
FRE  = 5. 

«*„«*««! 

C M=NUMRER  CE  RECCROS  IN  FILE 

C THR  E S = INTEGER  OATA  VALUE  USED  AS  THRESHOLD 

C SB=SPECTRAL  BAND  OF  INTEREST  l,2,CR  3 

C E I LE  = F I L E PROM  WHICH  OATA  IS  TAKEN  USED  ONLY 


FOR  LABELING 

t *a« •**«««* 


T=THRES 
N = 646 

CALL  PLOTSID.O.OI 

CALL  TCNEIO.  .0.. I RAT9.-16) 

0#*»**#****B*4flAPKEN  7 BY  4 INCH  AREA  TO  READY  PLOTEQ 
XI  1 I =0 • 

Y I 1 ) =0  . 

X(2)=0. 

Y<  2)  =700. 

XI 31=400. 

YC2>=700. 

XI  4)  =ACO 

Y ( 41=0. 

CALL  TONE  <X,Y,4,1) 

CALL  PICT i 0. .0. t-ocR) 
f*  = *?«  = *w  = = «*==LA3FL  RL  GTAS*** 

VALUE=T 

CALL  NIJmpER  (0.  ,0.  . 3*5 . , VALUE.0..-1  ) 

VAL  UE  = c 3 

CALL  NUMBER (0  . .50. , 3 5 . , VALUE . 0. , - 1 > 

CALL  NUMBER (0. ,lC0..35..FILE.0.,-t> 

CALL  °LCT  ( 100.  ,200.  . - 3 t 


<* 5~  ■?  = * e , 


!•>**«****»**»«**»■*»**** 


B = M 

CALL  PLOT  t B ,000. ,2) 

CALL  PLOTIB  ,646., 2) 

CALL  PLCTIO. ,646. , 2) 

CALL  PLOT ( 0. , 000. , 2) 

CALL  TONE ( 0. ,0. , IP AT9.-16 I 
►**»**» 

PRE  = 0 

DO  100  1=1 , M 
READ ( 2 ) ID 
C = 0 

00  10  A=SB, 1938,3 


rH1S  PaG£  IS  BEST  QUALITY  PRAUlIOARl  1 

"*-« <*« nMiau) io 
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?MID(  A t.LT.TI  GO  TO  5 
DATA ( C) “1 
GO  TO  10 
5 DATA ( C) »0 

10  CONTINUE 

DATA ( 64  71 »0 

IF(DATAUK£3.1>  GO  TO  20 
!F( PBE.EO.O)  GC  TO  100 
XI  1)«!-1. 

X( 2 » — I — 1 . 

XI  31 “I 
XI  41«! 

Y(1I«S 
Y( 2)*J-1. 

Yl 3)«J-1. 

Y(<tl*S 

CALL  TONE  (X.Y.A.ll 

oe.E*0 

GO  TO  100 

20  IF  (PBE.EQ.OI  00  Te  3C 

GO  TO  100 
30  5=J-l 

pp  = =1 

ioo  continue 

CALL  PL0T(0.t0.«*999l 

STOP 

ENO 


OUTPUT 


JOB  0900  9 V GOGGANS 


104  CAPOS  9E10  7 PAGES  POINTED 


This  pa  si  is  besi  quality  pbacii  mum 

too  M Ov>P  X yUK..  1 S HaL  j 0 DL  u — - 
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! 


c 


. 


C AREA  CONSIDERS  THE  SCANED  PICTURE  NATURE  OF  DIGITIZED 

C MULTISPECTR AL  SCANER  data  STORED  ON  MAGNETIC  TAPE.  AREA  GENERATES 

C A HISTOGRAM  OF  NUMBER  OF  PIXELS  ABGVC  A S ET  THRESHOLD  IN 

C CONTIGOUS  BLOCKS  OF  PIXELS  VS.  THEIR  RELATIVE  OCCURRENCE.  AREA 

C MAKES  THE  SAME  ASSUMPTIONS  ABOUT  THE  OATA  AS  PROGRAM  PLCTIT. 

C JOB  CONTROL  LANGUAGE  MUST  BE  INCLUDED  Tu  IDENTIFY  UNIT  2 AS 

C THE  TARE  CONTAINING  THE  FILE  TO  BE  USED. 

INTEGER  0ATAI646) , SB . GC ( l 000 ) . THRES 
DIMENSION  Xl( 1002) ,Yl( 10021 , 10(1938) 

INTEGER  P.O.P .S.T.U.X.Y.Z 
INTEGER  8(16000  ) ,CN( 160001 
DIMENSION  NGSI323I ,NGN< 32 3 ) . N3BI 323 ) 

INTEGER  OGS I 32 3 1 , OGN  t 323 ) .OGB I 323 ) 

C»—********t  PARAMETERS  WHICH  CHANGE  WITH  DATA  SETS  **•**•****»***••» 
THRES-184 
SB-3 
FILE-5. 

M-898 

C *»*•**•*•••«••*****•***«:••***••*••»—  ****•****•**•*•••••*••*•••••*••*< 

C THRES-INTEGER  DATA  VALUE  USED  AS  THRESHOLD 

C SB- SPECTRAL  BAND  OF  INTEREST  l, 2, OR  3 

C FILE-FILE  F RC m which  DATA  IS  TAKEN  USED  ONLY  FOR  LABELING 

C M-NUMSER  OF  RECORDS  IN  FILE 

C*«*****ft*  ******»•****»*»**—***»♦»*******••*•**—*•«**— ••*•**•*•••••« 

NC-0 

NB-0 

N-6A6 

DO  1 1-1,1000 

1 GC (I) -0 

00  60  I-UM 
RE AO ( 2)  ID 
NN-0 
K-l 

DO  3 J*  SB , 1938 , 3 
IF(ID( JI.LT.THRESI  GO  TO  2 
DATA ( K ) = l 
GO  TO  3 

2 - DAT  A( Kl -0 

3 K-K+l 
K-0 

00  8 J-I.N 

IF  (OATAl JI.EO.OI  GO  TO  7 

IF  (K.EO.l)  GO  TO  5 

NN-  NN»1 

NGSINNI-J 

NGN( NNI-l 

K*  1 

GO  TO  8 

5 NGN( NNl-NGNINNI+l 

GO  TO  8 

7 K=0 

8 CONTINUE 

IF  (NB.EO.OI  GO  TO  10 
00  9 Z-l.NB 

9 BtZJ»0 

10  IFINN.EO.OJ  GO  TO  35 

DO  11  R-l.NN 

11  NGBtRI-0 


«,sBESI»»“mm0n0ABU 
rmsPAGEis  ® 0BDC 
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DC  30  P*1.NN 
IF(NO.EO.O)  GO  TO  25 
DO  20  0* l v NO 

I F I NG  S( P I • GT. OGS I Q ) ) GO  TO  12 

U-0GSIQI-NGSIP1 

IF ( U.GT .NGN  I P ) ) GO  TO  20 

GO  TO  15 

U-NGSlP)-OGSm 

IFfU.GT.OGNIQI  I GO  TO  20 

IFINGBIP1.E0.0*  GO  TO  17 

IF  (NGB(PI.EO.OGBIO)I  GO  TO  20 

S-OGBIOI 

T-NGBIP) 

DO  16  R-l.NO 

IF ( CGBI  P.)  • NE. S ) GO  TC  16 

OGBtPI-T 

CONTINUE 

CNITI*  CN1T)*CN(SJ 
CNI S)«0 
GO  TO  20 
NGB ( P)=OGB<  0) 

S-OGBIO) 

CNI  S)«CNI S)*NGNIP» 

BI  S»»l 
CONTINUE 

IFINGB(PI.NE.O)  GO  TO  30 

NB«NR*l 

IF( NB.LT. 160001  GO  TO  27 
WRITEI6.26) 

FORMAT! *0', 'NUMBER  OF  BLOCKS  DIMENSION  EXCEDEO* » 
STOP 

NGB<P»-NB 

BINBI-l 

CNINBI-NGNIP) 

CONTINUE 

IFCNB.EO.O)  GO  TC  41 
DO  40  Y«1,NB 
IF(BIY).EO.l)  GO  TO  40 
IF  ICNIYI.EO.OI  GO  TO  40 
X-CNIY) 

1FfX.LT. 10001  GO  TO  37 
X-1000 

GCIXI-  GC I X I ♦ 1 
CNI Y)»0 
CONTINUE 

Irl NN.EO.O I GO  TO  55 
DC  50  Y-l.NN 
OGSI YI-NGSI V| 

OGNIYI-NGNIY) 

CG8I YI-NGBCYI 
NO-NN 
CONTINUE 

IF  INB.E0.01  60  TO  80 

DO  75  I-l.NB 
IFICNfl I.EO.OI  GO  TO  75 
Z-CNI I ) 

IFIZ.LT.1000  t GO  TO  70 
Z-1000 


fHlS  PAGE  IS  BEST  QUALITY  FRACIICAiUjl 
COPY  EUHi.ISHED  TO  QDS  , - 


93 


GC(Zt«GC<Z»«-l 

continue 

CONTINUE 
PCOUNT-O 
00  83  1-1,1000 
PCOUNT»PCOUNT«-GC(  I ) 

00  90  I-l, 1000 

xiin-i 

IFIGCII I.EO.O)  GO  TO  85 
VI ( I l-GC ( I l/PCOUNT 
GO  TC  «0 

vim-o 

CONTINUE 

NPS-200 

value-thres 

K-NPS+1 

DO  100  I -K  , 1000 

IF  I v 1 » I I.EO.O. I GO  TO  100 

WRITEI6.105I  I ,Yl(  I I,GC( I ) 

CONTINUE 

WRITEI6.U0)  PCOUNT 
WRITE(6,120)  S3, FILE 

FORMAT  ( • • * • I = Y!I)«  ',E9.2.'  Gill-  ',13* 

FORMAT! *0' , 'PCOUNT  « ',F9.2) 

FORMAT! 'O' ,' SB*  ',12,'  FILE-  '.FA.2) 

SB1-SB 

CALL  LI0NIXI,Y1,NPS,VALUE,F1LE,SB1) 

STOP 

ENO 


SUBROUTINE  LION! X, Y.NPS, VALUE, FILE, SB  I 

'•»»*»  — *LI  ON  GENERATES  ALL  AX  I S,  SC  Al  I NG,  LABELING  ANO  PLOT!'," 


NECESSARY  FOR  MAKING  HISTOGRAMS  -*=*»**»**»***»***- 
DIMENSION  X11002  l,Y!1002  1,1011938) 

CALL  PLOTSIO.O.O) 

CALL  PLOT ( 0. , 1 » ,2  I 

CALL  NUMBER  I 0. 0, 9. 5,. 35, VALUE, 0. 0,-1* 

CALL  NUMBER! 0. ,9.,.35,FILE,0.,-1> 

CALL  NUMREO(0.,8.5,.35»SB,0.,-1* 

CALL  PL0TI3., 1.25, -3* 

CALL  SCALE(Y,5.,NPS,l* 

CALL  SCALE!X,7.0,NRS,l) 

CALL  AXIS!0,,0., 'RELATIVE  CCCURR  ENC  E*  , «•!  9 , 5 .0,  90.  , Y ( NPS*  l I , 
3 Y( NPS»2I  ) 

CALL  AXISI0.,0., 'NUMBER  OF  PIXELS  IN  BLOCK', -25, 
a7.,00.,X(NPS+l I.XINPS+2*) 

CALL  LINEIX, Y.NPS, 1,0,0) 

CALL  PL0TH0..-1.25,-3I 
CALL  PLCT!0.,I.,2) 

CALL  PL0T(0.,0.,«-999) 

RETURN 

ENO 


OUTPUT 


JOB  8907  9V  GOGGANS 


188  CAROS  »EAD  8 PAGES  PRINTED 


0 CAROS  PUNCHED 


^HIS  PAGE  IS  BEST  QUALITY  PRACTICABL1 

CUl  ¥ EUmtlJHEjj  TO  DDG 


12 

15 


16 

17 


20 

25 


26  , 
27 


30 

35 


0C  30  P-l.NN 

IF  < NO. E 0*0  I GO  TO  25 

DO  20  0-1, NO 

IF1NGSIPI.GT.0GSIQII  GO  TO  12 
U-OGS(QI-NGSIP) 

IF ( U.GT .NGN! P ) ) GO  TO  20 

GO  TO  15 

U-NGSI P)— OGSI 01 

IFIU.GT.OGNIQl  ) GO  TO  20 

I F 1 NG8 ( P I . EQ. 0 I GO  TO  17 

IF  (NGBIPI.EO.OGB(OI)  GO  TO  20 

S-OG8IOI 

T-NGBIP1 

00  16  «-l, NO 

IF ( CGS<  PI. NE. S I GO  TC  16 

OGBIPI-T 

CONTINUE 

CN(T>»  CN I T 1 *CN! SI 
CNI S)»0 
GO  TO  20 
NGBI PJ-OGBI  0) 


S-OGBIOI 

t N ' ‘ ‘ ' 


:nisi«cn<si*ngnipi 

B(  SI-1 
CONTINUE 

IF  I NGBI PI .NE.O I GO  TO  30 
N3«NB*1 

IFINB.LT. 160001  GO  TO  2T 
MRI TE(6 ,26  I 

FORMAT! '0* , •NUMBER  OF  BLOCKS  DIMENSION  EXCEOEO* 1 
STOP 

NGB! PI-NB 
B(NSI*l 
CNINBI-NGNCPJ 
CONTINUE 

IF  I NB.EO. 01  GO  TC  61 
DO  60  Y-l.NB 
IFIBIVl.EC.il  GO  TO  60 
IF  ICNIYI.EO.OI  GO  TO  60 
X-CNIYI 

IF ( X.LT .10001  GO  TO  37 
X-1000 


37 

GC ! X I - GC ( X 1 + 1 
CNI Yl-0 

60 

CONTINUE 

61 

IF!  MN. EO.O 1 GO 
00  50  Y-l.NN 
OGSI YI-NGSI Y| 
0GNIY1-NGNIY1 

TO 

55 

50 

CGBI YI-NGB! Yl 

55 

NO-NN 

60 

CONTINUE 

IF  INB.EO.OI 

DO  75  I-l.NB 

IF  1 CNI I l.EO.OI 

GO 

TO 

80 

GO 

TO 

75 

Z-CNII) 

IFIZ.LT. 1000  1 
Z-1000 

GO 

TO 

70 

this  15 
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